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LIiprary 


[NESS AM SIN Gah 


A computer-aided Optimal Pilot Modeling Procedure was 
utilized to develop flight director laws for implementation 
in the longitudinal mode of a "head-down" cockpit display 
for a vertical take-off and landing aircraft. These flight 
director laws were intended to enhance the instrument hover 
and landing capabilities of the Bell UH-1H helicopter equipped 


with a multi-function stroke-written cathode ray tube display. 
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Pee LRODUCT ION 


Increasing emphasis is being placed on the development 
of vertical and short-field take-off and landing type air- 
craft. Of those that have been developed few, if any, have 
a satisfactory instrument capability. The conventional 
presentation of flight information is simply insufficient for 
the transition to hover flight, particularly in limited 
access landing situations. 

Current efforts to alleviate this deficiency have been 
most successful in the area of electronic displays. One 
Such display is the Hess V/STOLAND display developed at the 
Ames Research Center of the National Aeronautics and Space 
Administration by Professor R. A. Hess. This display pre- 
sents collocated conventional flight information for the 
initial phase of the approach and integrated contact situa- 
tion information for the transition and hover phase. Imple- 
mentation and evaluation by this author has resulted in 
various minor adjustments and modifications. For optimum 
effectiveness, however, a flight director system should be 
incorporated. A flight director system consists of display 
symbols which indicate desired Revenant CieEne DEloL' Ss 
controls, and the so-called flight director laws which drive 
these symbols as a function of aircraft deviations from a 


prescribed operating point. 





It is the purpose of this thesis to design the flight 
director laws for the longitudinal mode of the UH-1H 


helicopter in hover using the Hess V/STOLAND display. 
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tT, eee eeOPMENT OF THE HESS V/STOLAND DISPLAY 


It is the opinion of the author that a "head-down" VTOL 
landing display must provide a realistic representation of 
the actual visual world. This opinion was confirmed with 
the modification and evaluation of the Hess V/STOLAND 
display during an effort to enhance its capabilities in 
the hover and landing mode. 

The display was flown by the author in a fixed-based, 
digital simulation of the Bell UH-1H helicopter at the Ames 
Research Center. THe UH-lH is a single lifting rotor heli- 
copter with stabilizer bar and conventional tail rotor, but 
without a stability augmentation system. Hover and landing 
transitions from a straight-in approach were made in order 
to ascertain both the accuracy and the confidence of the 
pilot with the display. _ 

Presentations were made on a Sperry Flight System stroke- 
written multifunction display (MFD) with viSual references 
provided by TV simulation of an actual airport environment. 

Enclosure (5) of reference [2], an excerpt from 
"V/STOLAND MFD Display for Research Computer" by R. L. Sharp, 
represented the original display of Professor R. A. Hess. 
The hover mode of this display is shown in Figure 1. The 
ensuing development in [2] incorporated revisions by Pro- 
fessor Hess, and further changes made by the author during 


implementation. (Note that the flight director symbols 


Watt 





indicated were not operative.) This: display provided satis- 
factory results during the approach phase, but proved some- 
what inadequate during the transition and hover phase. A 
well deSigned flight director could alleviate this inadequacy. 
In the final analysis the Hess V/STOLAND display seemed 
highly adaptable for shipboard uSe due to its ability to 
present a dual perspective view of the landing situation, 
and through its potential to maintain the necessary flight 
control sensitivity for limited access landings through the 


use Of a flight director system. 
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Mtoe ve bOrMENT OF THE THEORETICAL MODEL 


A. PILOT MODELING 

The modeling hypothesis utilized for this development is 
that, subject to his inherent limitations, the well-trained 
well-motivated pilot behaves in an optimal manner. The 
pilot's control characteristics can be modeled by the solu- 
tion of an optimal linear control problem and optimal estima- 
tion problem with certain "modifications". These pilot 
modifications are: 

1. time delay - a pure time delay is included in each 
Siethe pilot's control outputs. 

2. neuromuscular dynamics = each output neuromuscular 
system is modeled as a first order lag. 

3. observation noise and motor noise - each observed 
variable is assumed to contain pilot induced additive white 
noise which scales with the variance of the observed variable; 
each control output is assumed to contain pilot induced addi- 
tive white noise which scales with the variance of the control 
Motion. 

4. if a variable is displayed explicitly, the pilot also 
perceives the first derivative of the variable but no higher 
derivatives. The first derivative is also noise contaminated. 

5. the index of performance for the optimization proce- 
dure iS chosen subjectively to mirror what the display system 
designer believes to be the task and control objectives as 


perceived by the pilot. 


ie) 





The time delay of the pilot is modeled through the use 
of a second order Pade' approximation, which permits direct 
implementation of the computer-aided optimal pilot modeling 
procedure. Table 1 presents the complete list of pilot 
modeling parameters. These parameter values, as well as 
the form of the observation and motor noise covariances have 


been well documented in the literature e.g. Refs. 1 and 3. 


Gee VEHICLE MODEL 

The vehicle model used was the unaugmented UH-1H heli- 
Peprer with a single lifting rotor with stabilization bar, 
a device attached to the rotor hub which provides pitch and 
Pomecdamnping, and a conventional tail rotor. The following 
assumptions were made: 

1. the vehicle was an idealized rigid airframe with a 
femerm attached. 

2. the rotor was described by the tip path plane whose 
orientation determines the propulsive and aerodynamic forces 
and moments. 

3. no rotor degrees of freedom are considered other 
than the control inputs which serve to describe the instan- 
taneous tip path plane orientation. 

4. all coupling between the longitudinal and lateral 
Meeron is ignored. 

5. linearized small perturbation motion about the 


horizontal reference flight path is assumed. 





The longitudinal equations of motion for the UH-1H are 


presented in Table 2. 


C. TURBULENCE MODEL 

immalreratt Stability and control analysis, gusts often 
play an important role. Accordingly, the most sophisticated 
and realistic treatment of gust disturbances should be em- 
ployed. Therefore, the stochastic approach was adopted in 
which the gust velocity in any direction at any fixed point 
in space is considered to be random in nature. Reference 
[4] suggest the use of the simplified spectrum for flight 
control system analysis presented in Table 3. There is no 
theoretical basis for this form, but its agreement with 
measured spectra, the preservation of the characteristic 
length, Lie and its simplicity suggest its use. For these 
reasons this spectrum was chosen for the development pre- 
sented here. It should be noted that for development of the 
hover mode flight director, only the horizontal turbulence 
will be considered. Strong rotor downwash, produced by high 
power requirements for hover flight in close ground proximity, 
tend to negate any effects of vertical turbulence. Also, 
Since the simplified spectrum form has been normalized to 
the reference airspeed (Un). hover flight would theoretically 
have no turbulence. This, in fact, is not the case; there- 
fore a reference airspeed of 20 knots (for turbulence modeling 


Only) has been introduced. 
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Dee PL LOT—-VEHIECLE SYSTEM MODEL 

The Pilot-Vehicle System Model depicted in Figure 2 
represents the optimal control and estimation model used 
for the development of this study, and Table 4 contains the 
system equations, As Figure 2 indicates, the vehicle is 
disturbed from the steady hover operating point by the hori- 
zontal turbulence. Linear functions of the system States 
are displayed to the pilot in the cockpit. The pilot then 
Operates upon these displayed variables and generates correc- 
tive control outputs. It is assumed that the manner in 
which the pilot operates on the displayed variables can be 
modeled by an optimal estimator and a state feedback con- 
troller. The optimal state feedback control, u(t), is then 
delayed, summed with the pilot white motor noise, and passed 
through the neuromuscular dynamics to form the pilot control 


meotLons. 
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Ev: 


titer lenis CONTROL—EoTIMATION PROBLEM 


Before discussing the specifics of the pilot modeling 


procedure, a review of the optimal control-estimation 


problem is in order. This problem can be outlined as 


follows: 


Consider a system described by 


x(t) = A x(t) + B u(t) + y w(t) 


pace) = C x(t) 
where 
A woman on on plant matrix 
fee tS an n x 1 state vector 
B is ann x p control matrix 
Mages ap x | control vector 
Y iS ann x t disturbance matrix 
w(t) is at x 1 disturbance vector 
meee is aq x 1 output vector 
g fomamces. 1 Output matrix 


PeGeme (t)eis @2 vector Of linearly uncorrelated, zero 


mean white noise Signals with Gaussian amplitude probability 


@stripution functions. The elements of w(t) are assumed to 


a7 





be sample functions from n random processes which are each 


ergodic and are jointly ergodic. The covariance matrix for 


w(t) is 
AL = 
pace) areata tT) ) = Fo6( 7) 


where 6(t) is the unit impulse function. 


The measured quantities Or sensor Signals are 


Gey meer te) ot) 
where 
(t) is a u x 1 meaSurement vector 
H Tease x nN Measurement matrix 


ie oeisuex Temeasumement nOlse vector 


moeselements of v(t) are assumed to be sample functions from 
p random processes each of which are ergodic and jointly 


eegodic. The covariance matrix for v(t) is 
a | 
aa even meet) G6 (+) 


The system is assumed to be completely controllable and 
completely observable. It is desired to find the control 
eem@ettOn H(t) which minimizes the quadratic scalar index 


of performance 
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eT ‘T 
Smee ow tye {t) sO y(t) + u(t) R u(t)) dt 


T'->00 0 
where 
Q is aq x q Symmetric output cost weighting matrix 
and at least positive semidefinite 
Peo canes xp symmetry ec Control cost weighting matrix 


and positive definite 


The solution to the linear quadratic Gaussian control problem 
can be outlined as follows: 

1. The optimization problem can, by the called Separa- 
tion Theorem, be broken up into two separate problems, an 
optimal control problem and an optimal estimation or filtering 
problem. 

2. The optimal estimation or filtering problem generates 
an optimal estimate, me) Speen otcaee em (t).)) Ihis estimate 
is optimal in the sense that 

At 


lim 1/T fs x(t) x(t) dt 
T+ 0 


1S minimized, where x(t) is the estimation error defined as 


x(t) = X(t) - x(t) 


The optimal estimator (or Kalman filter) has the form 
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The eStimator gains are given by 


K=PH G* 


mnere P Pomene error COVariance matrix 
= gel k 
ipo (t+ 1) )] = P 6(t) 


Meeece tne positive definite solution to the steady-state 


filter matrix Riccati equation 


[> 
[td 
+ 
[td 
Dy 
4+ 
[= 
‘es 
[x 
Has) 
Tt 
[Q 
ere 
[rd 
i 
| 


3. The optimal control problem generates an optimal 
SemerOl law u(t) which is a linear function of the estimated 


Scate 


ferme > iS ap x n Optimal controller gain matrix. The gain 
Mem@etx b> 1S identical to the one obtained by solving the 
Optimal control problem with no system disSturbance, exact 
State information, and the index of performance given by 

oO 


J Ss fy (t) Q@ y(t) + u(t) R ult)] at 
0 
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The controller gain matrix L is given by 


i, 2 mae 


where S is the positive definite solution to the steady- 


State control matrix Riccati equation 


It can be shown that the state covariance matrix 
Af 
Epesaey ex (€ + t)] = (P + M) O( 74) 


where P is the solution to the filter matrix Riccati equation 


Peet 2S the positive definite solution to 


In addition to the solutions outlined above, it can be 
shown that the transfer matrix relating the Laplace trans- 
form Of the Optimal control law u(t) to the Laplace transform 


of the measurement vector z(t) (with v(t) = 0) is given by 


(Ss) 


IN 


Eee ae eet Sat BLT KE) K 


2a 





u(s) =2 [u(t)] 


OS) Selene 


Ze 








Vee oon -VEntTCLE MODEL 


The equations of Table 4 are specialized as follows: 
the 13 differential equations describing the “plant" consist 
of the six vehicle state equations involving the state 


variables 


u - groundspeed 


Ww - velocity in the z stability axis 
arseer lon 

Gee) pitch rate 

6 -=- pitch angle 

h - altitude deviation 

x - longitudinal displacement from the pad 


In addition, a state We is added which represents the 
Mongitudinal turbulence. This Us represents the output of a 
shaping filter excited by white noise and has the power 
Spectral density given in Table 3. 

Four additional states A, > Cys were necessary for the 
second order Pade' approximation for the pilot's effective 
time delay for each control. 


Finally, two states 6, and 6 representing the pilot's 


(ead 
cyclic and collective control motions, respectively, were 
formed. These states resulted when the first order neuro- 


iiecular dynamics were put in state format. The A, B, and 


Y matrices are shown in Tables 5 and 6. 
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The Output variables, y(t) are 


ae =U 
7 
Y, = 

Wie aes 


itiete 7 Gives the corresponding C matrix. The displayed 
Variables, Zz, (t) are given in Table 8, while Table 9 gives 
Gite H matrix. 

Mee scOVaritance Of F and G are given in Table 10. Note 
Mitesall but one of the elements of F and G depend upon 
the variances of system variables. These variances, however, 
are not known a priori. Hence an iterative solution is 
necesSary. Initial guesses are made for the variances, a 
SOlution is obtained, the resulting variances used, and a 
new solution obtainec, etc. AS pointed out in Ref. [1], the 
number of iterations necessary for a final solution is 
roughly twice the number of displayed or perceived variables, 
in this case 2 x 6 = 12. 

meso ana k matrices are shown in Table 11. The elements 
of these matrices are the squares of the reciprocals of the 
“maximum allowable" deviations of the output quantities. 
Thus, when any output variable attains the allowable limits, 
it makes a contribution of unity to the integrand of the index 
Of performance. Table 12 lists these maximum allowable limits 
chosen subjectively for the particular vehicle and task under 


study. 
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The optimal control and estimation problem was solved 
Hsing a modified form of VASP (Ref. [5]) on the Naval 


Postgraduate School IBM 360 computer. 
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Vinee ool. OLRECTOR LAWS 


From the output of the final iteration of the modeling 
procedure, 12 pilot "transfer functions" were obtained. 
These 12 functions relate the Laplace transform of the cyclic 
and collective control motions to the transforms of the six 
displayed or perceived variables. 

The amplitude ratio plot of one Such transfer function 


6 
—=(s) 


is shown in Figure 3. Each of the 12 transfer functions 
exhibited similar frequency characteristics. The similarity 
included break frequencies and high frequency asymptotes. 
Pure gain approximations were made to the Bode diagrams as 
shown in Figure 3. These gains were normalized with respect 
to the largest magnitude in each set of six to obtain the 
flight director gains of Table 13. These gains can be used 


to generate the director laws according to 


a “+ ~ 
- Ky x K, Ww + BE. + K, q K 


yee ey 
x u gq B 


Bh h 


Nee Te Kee ek. ag + K. h + Kh 
¢ ey om Cg er Ch Cy 


The normalized Bode gains of the pilot transfer functions 


(Table 13) and the RMS values of the six vehicle motion 


26 





variables were then used to qualitatively assess the extent 
to which each motion variable affects each of the two pilot 
control outputs. Figures 4 and 5 were obtained by multi- 
plying the RMS values of each of the six vehicle motion 
variables displayed or perceived by the pilot by the Bode 
gains of the respective transfer function for the cyclic or 
collective control, and normalizing to the largest value 
obtained for each control. The bars in Figures 4 and 5 can 
be thought of as representing approximations to the relative 
amount of power in the cyclic or collective output which 
can be associated with the vehicle motion variable indicated. 
Figures 4 and 5 indicate that height deviation, pitch 
rate, height deviation rate, and pitch deviations, in order 
of importance are associated with cyclic control; while 
Prcen Late, height deviation rate, pitch, height deviation, 
and groundspeed deviations, in order of importance, are 


associated with collective control. 


ad 





Wile ae OuCiUoLONS 


Intuitively, one would assume that cyclic motion in a 
hover would primarily be determined by groundspeed, range, 
Pitch, and pitch rate in that order, with little effect 
shown by height deviation or height deviation rate. Figure 
4 shows however, that for this study, height deviation is 
of primary importance, with pitch rate, height deviation 
rate, and pitch making minor contributions. Similarly, one 
would expect collective motion to be determined by height 
deviation and height deviation rate primarily, but Figure 5 
presents pitch rate as being the primary Stimulus, with height 
deviation rate, pitch, and height deviation making major 
contributions, and groundspeed having a minor effect. 

Although the analytical results seem in opposition to 
those obtained by intuitive reasoning, only pilot-in-the-loop 
Simulation will provide definitive answers regarding the 


efficiency of the director laws. 
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VIII. RECOMMENDATIONS 


The flight director laws resulting from this study 
should be displayed under simulation and evaluated as to 
their effectiveness. For comparison, additional studies 
should be made uSing varied modeling parameters and weighting 
functions. One such set of weighting functions could be 
obtained by setting the maximum allowable values in Q and 
imecqual £O the respective RMS values obtained in the 


modeling procedure. 
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TABLE 1. 


Time delay 


Neuromuscular 
system 


Observation 
noise 


Motor noise 


Index of 
performance 


AU PEN DI xX A. 


TABLES 


Pilot Model Parameters 
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tT = 0.3 Secs. 
Ty = 0.2 secs. 
7 2 
Ve moma |Z. {t) } 
i Pi 
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2 Z an 
a, = p'mEfu;_,(t)]  (1=2,3) 
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di 
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TABLE 3 


Atmospheric Turbulence Model 


2 
= s oe Ns 1 
I, 2 
u 5 We 
Oe —z 
i 
u 
iL u 
u 
_ Ps 
L 
u 
= 600 ft. 


20 knots (for turbulence modeling only) 


Set sec, 
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TABLE 4 


System Equations for the Optimal Control 
and Estimation Problem 


fe) = A x(t) + B u(t) + y w(t) 
bare) = C x(t) 
eee = H x(t) + v(t) = 2, (t) ats) 


Elw(t) w°(ttt)] = F 6(t) 
T 
E[v(t) v (ttt)] = G 6(T) 
aera T 
Meer = 2 iy (t) O y(t) + u (t) Rult)] dt 
T+0 0 
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y Matrix 


al 
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C Matrix 
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